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Abstract.  We  present  the  results  of  an  intensive  coordinated  monitoring  campaign  in  the  optical  and  X-ray  wavelength  ranges  of 
the  low-mass,  pre-main  sequence  star  V410  Tau  carried  out  in  November  2001.  The  aim  of  this  project  was  to  study  the  relation 
between  various  indicators  for  magnetic  activity  that  probe  different  emitting  regions  and  would  allow  us  to  obtain  clues  on 
the  interplay  of  the  different  atmospheric  layers:  optical  photometric  star  spot  (rotation)  cycle,  chromospheric  Ho-  emission, 
and  coronal  X-rays.  Our  optical  photometric  monitoring  has  allowed  us  to  measure  the  time  of  the  minimum  of  the  lightcurve 
with  high  precision.  Joining  the  result  with  previous  data  we  provide  a  new  estimate  for  the  dominant  periodicity  of  V410Tau 
(1.87 1970  ±0.000010  d).  This  updated  value  removes  systematic  offsets  of  the  time  of  minimum  observed  in  data  taken  over  the 
last  decade.  The  recurrence  of  the  minimum  in  the  optical  lightcurve  over  such  a  long  timescale  emphasizes  the  extraordinary 
stability  of  the  largest  spot.  This  is  confirmed  by  radial  velocity  measurements:  data  from  1993  and  2001  fit  almost  exactly 
onto  each  other  when  folded  with  the  new  period.  The  combination  of  the  new  data  front  November  2001  with  published 
measurements  taken  during  the  last  decade  allows  us  to  examine  long-term  changes  in  the  mean  light  level  of  the  photometry 
of  V410  Tau.  A  variation  on  the  timescale  of  5.4  yr  is  suggested.  Assuming  that  this  behavior  is  truly  cyclic  V410  Tau  is  the  first 
pre-main  sequence  star  on  which  an  activity  cycle  is  detected.  Two  X-ray  pointings  were  carried  out  with  the  Chandra  satellite 
simultaneously  with  the  optical  observations,  and  centered  near  the  maximum  and  minimum  levels  of  the  optical  lightcurve. 

A  relation  of  their  different  count  levels  to  the  rotation  period  of  the  dominating  spot  is  not  confirmed  by  a  third  Chandra 
observation  carried  out  some  months  later,  during  another  minimum  of  the  1.87  d  cycle.  Similarly  we  find  no  indications  for 
a  correlation  of  the  H<r  emission  with  the  spots’  rotational  phase.  The  lack  of  detected  rotational  modulation  in  two  important 
activity  diagnostics  seems  to  argue  against  a  direct  association  of  chromospheric  and  coronal  emission  with  the  spot  distribution. 

Key  words,  stars:  individual:  V410Tau  -  stars:  late-type,  coronae,  activity  -  X-rays:  stars 


1.  Introduction 

V410Tau  is  an  analog  for  the  young  Sun  on  the  pre-main  se¬ 
quence  (PMS).  Due  to  the  lack  of  strong  emission  lines  and 
infrared  excess  it  can  be  classified  as  a  weak-line  T  Tauri  star 
(wTTS).  This  term  defines  PMS  stars  without  obvious  signs  for 
disk  accretion,  while  young  stars,  in  which  accretion  from 
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a  circumstellar  disk  is  responsible  for  ultraviolet  and  infrared 
excess  emission  and  for  a  moderate  to  strong  emission  line 
spectrum  superimposed  on  the  photospheric  spectrum,  are 
called  classical  T  Tauri  stars  (cTTS).  Indeed,  V410Tau  shows 
a  small  infrared  excess  which,  however,  is  attributed  to  one  or 
two  close  companions  at  sub-arcsecond  separation  (Ghez  et  al. 
1993;  Ghez  et  al.  1997). 

WTTS  represent  an  evolutionary  stage  where  the  disk  has 
already  dissipated,  and  therefore  their  variability  is  not 
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related  to  accretion  processes  but  believed  to  be  a  manifesta¬ 
tion  of  magnetic  activity,  similar  to  that  observed  in  the  Sun, 
but  enhanced  by  several  orders  of  magnitudes.  Magnetic  activ¬ 
ity  seems  to  explain  all  the  optical  variability  phenomena  dis¬ 
covered  on  V410Tau,  comprising  a  large  range  of  time  scales: 
i )  variability  of  the  amplitude  of  the  optical  light  curve  on 
time  scales  of  years;  ii)  changes  in  the  Hu  profile  and  intensity 
recorded  within  months;  iii)  repeated  modulation  of  the  optical 
brightness  with  a  period  of  1.87  d;  and  iv)  sudden  brightenings 
and/or  emission  line  enhancements  evolving  on  time  scales  of 
hours. 

Items  i)  and  iii)  point  directly  at  the  presence  of  star 
spots.  The  periodic  1.87d  variability  of  the  optical  lightcurve 
of  V410Tau  was  first  reported  by  Rydgren  &  Vrba  (1983), 
and  confirmed  by  subsequent  observations  (Vrba  et  al.  1988; 
Bouvier  &  Bertout  1989).  It  is  attributed  to  stellar  spots  which 
produce  a  periodic  rotation  pattern.  From  these  studies  it  is 
shown  that  about  15-45%  of  the  stellar  surface  of  V410Tau 
should  be  covered  by  spots  in  order  to  explain  the  large  am¬ 
plitudes  observed  in  the  optical  bands.  The  models  indicated 
temperatures  for  the  spot(s)  that  are  500  K  to  1200  K  cooler 
than  the  photosphere,  similar  to  the  solar  case  (Wallace  et  al. 
1996).  Independent  evidence  confirming  the  spot  hypothesis 
was  provided  by  high-resolution  spectroscopic  observations: 
large  (>10%)  fractions  of  the  stellar  surface  at  about  1000K 
below  the  photospheric  temperature  can  considerably  alter  the 
profile  of  photospheric  lines  (Vogt  &  Penrod  1983;  Vogt  et  al. 
1987).  The  first  Doppler  images  of  V410Tau  were  published 
by  loncouret  al.  (1994)  and  Strassmeier  et  al.  (1994),  and  were 
confirmed  by  Hatzes  (1995)  and  Rice  &  Strassmeier  (1996). 
All  of  them  show  a  large  cool  spot  that  extends  over  one  of  the 
stellar  poles  and  smaller  cool  features  located  at  low  latitudes. 

Longterm  variability  of  the  photometric  amplitude  could  be 
related  to  changes  in  the  number,  size  and  location  of  spots. 
Longterm  periodic  changes  in  photometric  lightcurves  of  a 
number  of  solar-analogs  (Messina  &  Guinan  2002;  Berdyugina 
et  al.  2002)  and  more  evolved  RSCVn  binaries  (e.g.  Henry 
et  al.  1995;  Rodono  et  al.  2000;  Olah  et  al.  2000)  have 
been  interpreted  as  magnetic  activity  cycles.  Among  younger 
PMS  stars  a  lack  of  dedicated  monitoring  programs  has  long 
impeded  any  systematic  investigation  of  this  issue.  First  results 
on  longterm  photometry  including  PMS  stars  were  reported  by 
Strassmeier  et  al.  (1997)  based  on  observations  at  automated 
photoelectric  telescopes.  Thanks  to  the  efforts  of  Vrba  et  al. 
(1988),  Herbst  (1989),  Petrov  et  al.  (1994),  Strassmeier  et  al. 
(1997),  and  Grankin  (1999)  V410  Tau  has  by  far  the  largest  data 
base  of  photometry  among  the  PMS  stars.  Although  smooth 
changes  in  the  amplitude  of  the  optical  light  curve  have  been 
reported,  no  cyclic  or  predictable  behavior  has  been  found  yet. 

The  origin  of  variations  in  emission  lines  [items  ii)  and  iv) 
of  the  above  list]  is  sometimes  less  obvious.  Ha  observations 
of  V410  Tau  reveal  slow  changes  both  in  the  emission  line  pro¬ 
file  and  its  intensity.  There  are  two  kinds  of  changes:  smooth 
variations  of  the  line  profile  and  intensity  within  weeks  (e.g. 
Petrov  et  al.  1994),  that  sometimes  correlate  with  the  photo¬ 
metric  1.87  d  period  (Fernandez  &  Miranda  1998)  and  other, 
stronger,  variations  that  are  noticed  when  comparing  data 
taken  months  apart  from  each  other  (previous  references; 


Hatzes  1995).  Petrov  et  al.  (1994)  suggest  a  chromospheric 
origin  for  the  narrow  emission  peak  that  is  often  superim¬ 
posed  on  a  wide  platform  (see  also  Hatzes  1995;  Fernandez 
&  Miranda  1998),  while  the  platform  or  flat  extended  emission 
wings  might  arise  from  the  circumstellar  gas  environment. 

The  different  line  fluxes  between  observations  taken  at  dif¬ 
ferent  epochs  may  be  due  to  flaring.  Flares  are  magnetic  re¬ 
connection  events,  and  can  be  as  short  as  ~1  h,  such  that  their 
evolution  is  difficult  to  trace  in  non-continuous  observations. 
For  this  reason  the  number  of  flares  reported  for  V410Tau  is 
quite  limited.  Most  of  these  events  were  observed  spectroscop¬ 
ically  (e.g.  Welty  &  Ramsey  1995),  and  some  flares  were  de¬ 
tected  also  in  photometric  monitoring  programs  (Rydgren  & 
Vrba  1983;  Vrba  et  al.  1988). 

In  addition  to  these  well-established  activity  phenomena 
observed  in  the  optical,  magnetic  processes  further  up  in  the 
atmosphere,  i.e.  the  corona,  are  expected  to  give  rise  to  radio 
and  X-ray  emission.  Indeed,  the  radio  emission  of  V410Tau 
is  highly  variable  both  in  intensity  and  spectral  index  (Cohen 
et  al.  1982;  Bieging  et  al.  1984),  leading  to  the  conclusion  that 
it  must  arise  from  a  non-thermal  process,  probably  related  to 
magnetic  activity  (see  Bieging  &  Cohen  1989  and  references 
therein).  Bieging  &  Cohen  (1989)  monitored  the  radio  flux  den¬ 
sity  at  monthly  intervals  over  one  year.  They  failed  to  detect 
strong  radio  flares,  but  found  evidence  for  a  modulation  with 
a  period  that  is  half  of  the  period  reported  from  optical  ob¬ 
servations.  X-ray  data  of  V410Tau  has  been  presented  by  e.g. 
Strom  &  Strom  (1994),  Neuhauser  et  al.  (1995),  and  Stelzer 
&  Neuhauser  (2001).  Costa  et  al.  (2000)  has  extensively  dis¬ 
cussed  archived  International  Ultra-violet  Explorer  (IUE)  and 
ROSAT  observations  of  this  star.  Despite  its  obvious  variability 
in  the  X-ray  regime  no  direct  signs  for  rotational  modulation  or 
X-ray  flares  have  been  observed  from  V410Tau  so  far. 

The  detailed  relations  among  the  various  atmospheric 
regions  involved  in  stellar  magnetic  activity  (photosphere, 
chromosphere  and  corona)  remain  unexplained  and  call  for  a 
systematic  investigation.  Towards  this  end  we  organized  con¬ 
temporaneous  optical  and  X-ray  observations  of  V410Tau. 
Dedicated  photometric  as  well  as  intermediate-  and  high- 
resolution  spectroscopic  observations  in  the  optical  were 
planned  for  a  time  interval  of  ~  1 1  d,  simultaneous  to  three  ex¬ 
posures  with  the  Chandra  X-ray  satellite.  This  enabled  us  to 
study  activity  on  V410Tau  that  evolves  on  short  timescales, 
i.e.  rotational  effects  and  flares.  In  addition  we  examined  his¬ 
torical  photometric  data  of  V410Tau  to  examine  its  long-term 
variability. 

We  will  present  our  results  in  a  series  of  two  papers.  In 
the  present  paper  we  focus  on  the  examination  of  the  relation 
between  emission  in  various  energy  bands  and  the  optical  rota¬ 
tion  cycle,  and  the  long-term  photometry,  while  in  a  subsequent 
paper  (Fernandez  et  al.,  in  prep.)  we  will  concentrate  on  the  nu¬ 
merous  flares  found  during  our  campaign.  The  present  paper  is 
structured  as  follows.  In  Sect.  2  we  present  the  layout  of  our  ob¬ 
servations.  We  discuss  the  evolution  of  the  photometric  1.87  d 
cycle  over  the  years  in  Sect.  3.  This  section  includes  a  new 
determination  of  this  period  and  the  ephemeris  for  the  mini¬ 
mum  brightness,  and  the  tentative  detection  of  an  activity  cy¬ 
cle.  The  relation  between  spots  and  various  activity  diagnostics 
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Fig.  1.  Observing  Journal  for  our  monitoring  of  V410Tau  in 
Nov.  2001.  Note,  that  our  campaign  includes  a  third  Chandra  obser¬ 
vation  which  was  performed  in  March  2002  and  is  not  shown  in  this 
diagram. 


are  examined  in  Sect.  4.  We  devote  Sect.  5  to  the  description  of 
the  X-ray  spectrum  of  V410Tau.  The  results  are  discussed  in 
Sect.  6.  Section  7  presents  a  brief  summary. 


2.  Observations  and  data  reduction 

Simultaneous  optical  and  X-ray  observations  were  planned 
for  V410Tau  from  Nov.  15  to  26,  2001  (UT).  Optical  pho¬ 
tometry  and  spectroscopy  were  carried  out  from  several  ob¬ 
servatories  (see  below)  while  three  X-ray  observations  were 
scheduled  with  Chandra  using  the  Advanced  CCD  Imaging 
Spectrometer  for  Spectroscopy  (ACIS-S).  Our  observing  strat¬ 
egy  was  to  maximize  the  probability  of  detecting  large  ampli¬ 
tude  variability  in  the  X-ray  observations.  If  the  X-ray  emission 
comes  from  spotted  regions  the  maximum  of  the  X-ray  emis¬ 
sion  should  occur  at  the  minimum  of  the  optical  lightcurve,  i.e. 
at  phases  when  the  spot  is  on  the  visible  hemisphere,  and  vice 
versa.  Therefore,  the  Chandra  observations  were  scheduled  for 
the  time  near  minimum  and  maximum  of  the  known  rotation 
cycle,  respectively. 

Nevertheless,  the  third  Chandra  observation  was  de¬ 
layed  by  several  months  because  at  the  anticipated  time  in 
November  2001  the  satellite  had  to  be  shut  down  due  to  high 
solar  activity.  Therefore,  optical  data  is  available  only  for  the 
first  two  Chandra  observations.  The  complete  observing  log 
for  the  campaign  in  Nov.  2001  is  given  in  Fig.  1.  In  the  fol¬ 
lowing  the  individual  observations  and  the  data  reduction  are 
described. 


2.1.  Optical  photometry 

2.1 .1 .  Sierra  Nevada 

The  photometric  observations  carried  out  at  the  Observatory 
of  Sierra  Nevada  (Granada,  Spain)  were  performed  with  a 
Stromgren  photometer  attached  to  the  90  cm  telescope.  The 
photometer  is  equipped  with  identical  six-channel  uvby  spec¬ 
trograph  photometers  for  simultaneous  measurements  in  uvby 
or  the  narrow  and  wide  H/3  channels  (Nielsen  1983).  We 
carried  out  uvby  measurements.  Integration  times  were  60s 
for  V410Tau  and  for  its  two  comparison  stars,  HD  27159  and 
HD  283561,  as  well  as  for  the  background  sky.  Differential 
photometry  was  carried  out  during  all  the  photometric  nights, 
except  during  the  last  one,  in  which  the  absolute  calibration 
was  done. 

The  error  bars  of  the  differential  photometry,  as  esti¬ 
mated  from  the  difference  between  the  two  comparison  stars, 
are  0.005  mag  in  the  u  band  and  less  than  0.001  mag  for 
the  vby  bands.  Nevertheless,  we  note  that  both  comparison 
stars  are  about  two  magnitudes  brighter  than  V410Tau,  there¬ 
fore,  slightly  larger  errors  are  expected  for  the  target  star.  For 
the  absolute  photometry,  the  residuals  of  the  standard  stars 
(difference  between  computed  and  published  values)  are  less 
than  0.026  mag  for  the  V  band  and  less  than  0.008,  0.010 
and  0.017  mag  for  the  ( b  -  if),  ml  and  cl  indexes,  respectively. 

2.1 .2.  Mt.  Maidanak 

Diaphragm  photometry  was  performed  at  the  Mt.  Maidanak 
Observatory,  Uzbekistan.  All  observations  were  obtained  with 
the  48  cm  telescope  equipped  with  a  one-channel  pulse¬ 
counting  photometer.  We  carried  out  Johnson  UBVR  measure¬ 
ments  with  a  28"  diaphragm.  A  detailed  description  of  the 
equipment  and  measuring  techniques  was  given  by  Shevchenko 
(1980).  The  exposure  time  was  typically  40- 120  s,  depending 
on  the  band  and  sky-background  brightness.  As  a  rule,  seven 
standard  stars  from  the  list  of  Landolt  (1992)  were  observed 
on  each  night  to  determine  the  atmospheric  extinction  coeffi¬ 
cients.  The  instrumental  photometric  system  was  reduced  to 
the  standard  system  by  the  method  of  Nikonov  (1976).  The  star 
BD+28  643  was  used  as  a  secondary  standard.  The  rms  error 
of  a  single  measurement  on  a  moonless  night  was  0.032, 0.007, 
0.005,  and  0.005  mag  in  U ,  B,  V,  and  R,  respectively. 

2.1.3.  USNO  Flagstaff 

Observations  from  the  USNO,  Flagstaff  Station  (NOFS,  USA) 
used  the  NOFS  1.0  m  f/7.3  R/C  telescope  along  with  a 
SITe/Tektronix  2048x2048  CCD.  A  mask  around  the  CCD  lim¬ 
its  its  usable  area  to  1815xl815pixels,  or  about  20x20'.  By  po¬ 
sitioning  V410Tau  in  one  corner  of  the  array  the  largest  num¬ 
ber  of  candidate  comparison  stars  could  be  included.  Several 
nights  prior  to  the  campaign,  plus  the  photometric  nights  dur¬ 
ing  the  campaign,  were  used  to  calibrate  all  potential  compar¬ 
ison  stars  in  the  field.  UBVRCIC  filters  were  used,  along  with  a 
large  number  of  Landolt  standard  stars  (see  Landolt  1983  and 
Landolt  1992),  to  calibrate  the  comparison  stars.  Five  of  the 
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potential  comparison  stars  were  found  to  be  constant  during  the 
survey  interval,  and  bright  enough  to  give  reasonable  signal-to- 
noise  (S/N)  for  differential  photometry.  Differential  photome¬ 
try  at  UBVRCIC  was  performed  extensively  on  the  five  nights 
Nov.  16  through  Nov.  20  UT,  plus  a  few  data  points  on  nights 
prior  to  and  subsequent  to  this  period. 

2.2.  Optical  spectroscopy 

2.2.1.  Intermediate-resolution  spectroscopy 

Intermediate  resolution  spectroscopy  was  carried  out  at  the 
1.5  m  telescope  on  the  Observatory  of  Sierra  Nevada  (Granada, 
Spain)  using  the  spectrograph  ALBIREO  (see  Sanchez  et  al. 
2000)  in  long-slit  mode.  The  detector  was  a  EEV8821/T  CCD 
of  1152  x  770  pixels  with  a  22.5  pm  pixel  size.  Observations 
were  done  on  two  spectral  ranges:  4000-5160  A  (blue) 
and  5645-6790  A  (red).  The  full  width  half  maximum 
(FWHM)  of  the  lines  of  the  calibration  lamps  were  1.7  A 
and  1.5  A  for  the  blue  and  red  ranges,  respectively.  Several 
spectrophotometric  standard  stars  were  observed  on  the  three 
photometric  nights,  with  the  aim  of  correcting  the  slope  of 
the  continuum  but  no  flux  calibration.  A  set  of  templates,  with 
spectral  types  ranging  from  K2  to  K5  were  also  observed. 

The  data  reduction  and  analysis  was  done  with  the  longslit 
package  of  IRAF  (Image  Reduction  and  Analysis  Facility1). 
The  longslit  tasks  were  required  due  to  the  strong  distortion  of 
the  sky  lines  along  the  spatial  direction.  Individual  spectra  have 
a  S/N  of  about  30. 

2.2.2.  High-resolution  spectroscopy 

High-resolution  spectroscopic  observations  were  performed 
using  the  Fiber  Optics  Cassegrain  Echelle  Spectrograph 
(FOCES)  and  a  2048  x  2048  pixel  SITe  CCD  attached  to  the 
2.2  m  telescope  at  the  Calar  Alto  Observatory  (Almeria,  Spain). 
The  spectra  include  some  seventy  orders,  covering  a  range 
from  4200  A  to  7000  A,  with  a  nominal  resolving  power  of 
T/AT  *  30000.  Technical  details  regarding  FOCES  can  be 
found  in  Pfeiffer  et  al.  (1998). 

Additionally,  high-resolution  spectra  were  taken  at  the 
University  of  California’s  Tick  Observatory  on  Mt.  Hamilton 
using  the  Hamilton  echelle  spectrograph  on  the  3m-Coude 
telescope.  The  instrument  yielded  107  spectral  orders  span¬ 
ning  a  wavelength  range  of  -3500-10000  A.  A  thinned  Ford 
2048  X  2048  pixel  CCD  with  15  pm  pixel  size  was  used. 
Using  an  aperture  plate  above  the  slit  gave  a  slit  width 
of  640  pm,  corresponding  to  1.2"  projected  on  the  sky, 
and  ~2  pixels  on  the  CCD.  This  gave  a  2-pixel  spectral  reso¬ 
lution  of  -0.1  A  FWHM  (i.e.,  2-pixel  R*  60  000). 

All  high-resolution  spectra  have  been  reduced  and  ex¬ 
tracted  using  the  standard  IRAF  reduction  procedures  (bias 
subtraction,  flat-field  division  and  optimal  extraction  of  the 


1  IRAF  is  distributed  by  the  National  Optical  Astronomy 
Observatories,  which  is  operated  by  the  Association  of  Universities 
for  Research  in  Astronomy,  Inc.  (AURA)  under  cooperative  agree¬ 
ment  with  the  National  Science  Foundation. 


Table  1.  Observing  log  for  the  three  Chandra  ACIS-S  observations  of 
V410Tau.  Rotational  phases  are  computed  using  the  new  ephemeris 
given  in  Table  2.  Values  refer  to  the  0.2-8  keV  energy  band. 


Seq.# 

Start  Date 

[d/m/y  h:m] 

Exp.  Time 

[s] 

0rot 

Rate 

[10~3  cps] 

200130 

16/11/01  20:11 

14925.3 

0.91-1.00 

296  ±4 

200190 

19/11/01  23:37 

11007.6 

0.58-0.65 

361  ±6 

200191 

07/03/02  06:16 

17734.2 

0.89-0.99 

324  ±4 

spectra).  Wavelength  calibration  was  done  using  spectra  of  a 
Th-Ar  lamp.  Finally,  for  each  order  we  normalized  the  ex¬ 
tracted  spectra  by  means  of  a  polynomial  fit  to  the  observed 
continuum. 

2.3.  X-rays 

The  Chandra  observations  were  performed  with  V410Tau  on 
the  back-illuminated  ACIS-S3  CCD.  In  order  to  avoid  pile-up 
(the  detection  of  more  than  one  photon  as  a  single  event  which 
leads  to  distortion  of  the  spectral  shape  and  underestimate  of 
the  count  rate)  the  source  was  placed  4'  off  the  aimpoint.  This 
way,  for  the  expected  X-ray  brightness  of  V410Tau  in  the 
Chandra  energy  band,  pile-up  should  be  below  10%.  However, 
we  performed  additional  checks  on  the  count  rate  and  the  spec¬ 
trum  to  verify  that  pile-up  is  negligible  (see  Sect.  5).  Our  data 
analysis  is  based  on  the  events  level  1  data  provided  by  the 
pipeline  processing  at  the  Chandra  X-ray  Center  (CXC),  and 
was  carried  out  using  the  CIAO  software  package2  version  2.3. 

In  the  process  of  converting  the  level  1  events  file  to  a 
level  2  events  file  for  each  of  the  observations  we  performed 
the  following  steps:  We  filtered  the  events  file  for  event  grades 
(retaining  the  standard  ASCA  grades  0,  2,  3,  4,  and  6),  and  ap¬ 
plied  the  standard  good  time  interval  (GTI)  file.  Events  flagged 
as  cosmic  ray  afterglow  were  retained  after  inspection  of  the 
images  revealed  that  a  substantial  number  of  source  photons 
erroneously  carry  this  flag.  We  also  checked  the  astrometry 
for  any  known  systematic  aspect  offsets  using  CIAO  soft¬ 
ware,  and  performed  the  necessary  corrections  by  updating  the 
FITS  headers. 

To  find  the  precise  X-ray  position  of  V410Tau  we  ran  the 
wavdetect  source  detection  routine  (Freeman  et  al.  2002).  For 
the  source  detection  we  used  a  binned  image  of  the  ACIS-S3 
chip  with  ~3"  X  3"  pixels.  The  significance  threshold  was  set 
to  1 0  6,  and  wavelet  scales  between  1  and  8  were  applied. 
After  V410  Tau  was  located  in  this  way  we  extracted  the  counts 
from  a  circular  region  of  4"  radius  around  the  position  found 
by  wavdetect.  This  area  includes  ~91%  of  the  source  photons 
(for  a  monochromatic  1.49  keV  source).  The  background  was 
extracted  from  an  annulus  in  a  source-free  region  surround¬ 
ing  V410Tau.  The  observing  log  for  the  three  Chandra  expo¬ 
sures  and  the  average  ACIS-S  count  rates  of  V410Tau  in  the 
0.2-8  keV  energy  range  are  given  in  Table  1. 

2  CIAO  is  made  available  by  the  CXC  and  can  be  downloaded  from 
htp : //cxc . harvard.edu/ciao/download-ciao-reg . html 
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Table  2.  Ephemeris  of  V410Tau:  value  by  P94  and  our  update. 


Reference 

Tmin  [JD] 

P rot  [d] 

P94 

2446  659.4389 

1.872095  ±  0.000022 

here 

2452  234.285971 

1.871970  ±0.000010 

3.  Variability  in  the  optical 

The  rotational  period  of  V4 10  Tau  was  first  determined  by 
Rydgren  &  Vrba  (1983),  who  found  a  value  of  1.92  d  from  pho¬ 
tometric  data  obtained  over  a  6-day  interval.  Vrba  et  al.  (1988) 
combined  data  taken  during  five  observing  seasons  between 
1981  and  1987,  and  improved  the  measurement  of  the  period. 
Only  one  year  later  Herbst  (1989)  noted  a  linear  trend  present 
in  the  O-C  (observed  minus  computed)  diagram  for  minimum 
light  of  V4 10  Tau  based  on  the  ephemeris  by  Vrba  et  al.  (1988) 
indicating  the  need  for  a  further  revision.  The  most  recent 
and  most  widely  used  ephemeris  for  the  photometric  mini¬ 
mum  of  V410Tau  is  the  one  presented  by  Petrov  et  al.  (1994): 
JD(min)  =  JD  2  446  659.4389  +  1.872095(±0.000022)  E.  This 
result  is  based  on  data  from  1986-1992. 

We  have  used  this  period  and  ephemeris  to  phase-fold  the 
photometry  from  November  2001,  and  found  that  the  minimum 
is  offset  from  phase  0  =  0  by  A0  ~  0.14  (Stelzer  et  al.  2002). 
A  shift  of  the  optical  minimum  with  respect  to  the  ephemeris 
given  by  P94  had  already  been  identified  by  Grankin  (1999) 
in  his  investigation  of  seasonal  lightcurves  acquired  within  the 
last  decade.  Our  measurement  continues  the  monotonic  trend 
observed  since  1990.  This  migration  of  the  minimum  could  ei¬ 
ther  indicate  a  change  in  the  latitude  of  the  spots  or  the  need 
for  a  new  estimate  of  the  period.  The  structure  of  the  lightcurve 
of  V410Tau  undergoes  changes  over  the  years  (Fig.  3  and 
Sect.  3.3)  suggesting  that  indeed  the  distribution  of  surface  fea¬ 
tures  is  variable  in  time.  However,  spot  models  based  only  on 
photometric  data  cannot  provide  unique  information  about  the 
shape  and  latitude  of  the  spots. 

3.1.  Migrating  spots  and  differential  rotation 

As  outlined  above,  the  systematic  phase  shift  of  the  mini¬ 
mum  observed  over  the  last  ~10yrs  could  be  due  to  a  lat¬ 
itudinal  migration  of  the  spots  on  the  differentially  rotating 
surface.  The  differential  rotation  of  V410Tau  is  known  from 
Doppler  Imaging,  where  the  differential  rotation  parameter 
k  —  AQ/ncq  (Qeq  is  the  rotation  rate  at  the  equator)  for  an 
assumed  solar-like  rotation  law  was  found  to  be  ~ 0.001  (Rice 
&  Strassmeier  1996).  This  value  is  much  smaller  than  the  so¬ 
lar  differential  rotation  which  is  =0.2  from  equator  to  pole 
(Snodgrass  1983),  but  not  untypical  for  a  fast  rotating  star  (see 
e.g.  Collier  Cameron  2002). 

The  phase  shift  accumulated  over  the  last  decade  with  re¬ 
spect  to  the  ephemeris  by  P94  is  A 0  =  0.14  (see  Stelzer  et  al. 
2002).  This  indicates  that  A  P/P  ~  6  x  1 0  "5.  If  we  assume  a 
solar  rotation  law,  and  that  the  major  spot  was  centered  on  the 
pole  initially  (0\  —  90°),  this  change  in  period  implies  that  the 
spot  should  have  moved  to  a  latitude  of  62  ~  75°  by  2001. 


However,  such  a  large  displacement  is  clearly  inconsistent  with 
the  relatively  unchanged  shape  of  the  photometric  lightcurve. 
Furthermore,  Doppler  maps  obtained  in  the  years  1990  to  1994 
seem  not  to  show  a  systematic  movement  of  the  main  spot. 

3.2.  An  update  of  the  period  and  ephemeris 

In  this  section  we  examine  the  second  possibility  for  the 
observed  systematic  shift  of  the  minimum  in  the  optical 
lightcurve,  namely  the  need  for  a  refinement  of  the  value  for 
the  rotational  period.  In  the  following,  whenever  we  refer  to  the 
“rotational  period"  or  “rotation  cycle”  the  reader  should  keep 
in  mind  that  this  is  the  period  representative  for  the  main  spot 
that  dominates  the  observed  lightcurve. 

For  our  update  of  the  period  we  made  use  of  historical 
data  from  the  T  Tauri  photometry  data  base3.  We  split  the 
available  V  band  photometry  on  a  yearly  basis,  and  derived 
the  most  significant  period  for  each  of  the  observing  intervals 
using  the  string  length  method  (Dworetsky  1983).  Each  sea¬ 
sonal  lightcurve  was  folded  with  its  period  to  determine  the 
time  of  minimum  To  by  fitting  a  polynomial  to  the  folded 
lightcurve.  We  used  the  value  for  the  period  given  by  P94  to 
compute  the  number  of  cycles  elapsed  between  each  of  the 
seasonal  points  To  and  the  time  of  minimum  observed  by  us 
in  November  2001,  using  this  last  observation  as  a  reference 
point.  Then  we  computed  the  O-C  diagram  for  the  minima  To 
as  a  function  of  cycle  number.  The  residuals  in  this  diagram 
can  be  minimized  by  modifying  the  period.  But  the  O-C  resid¬ 
uals  show  systematic  non-linear  trends  for  data  acquired  be¬ 
fore  1990.  The  large  scatter  for  these  earlier  years  may  indi¬ 
cate  that  the  spots  migrated  in  an  irregular  way  on  the  surface 
of  V410  Tau.  For  years  later  than  1990  we  observe  a  monotonic 
trend  in  the  O-C  residuals.  Therefore  we  used  only  observa¬ 
tions  obtained  in  1991  and  later  for  the  final  determination  of 
the  period.  This  way  a  best  fit  period  of  1.87 1970(10)  d  is  found 
(see  Table  2  for  the  new  ephemeris  of  V410  Tau),  slightly  lower 
than  the  earlier  determination  by  P94  which  was  based  on  data 
from  1986  to  1992.  The  residuals  in  the  O-C  plot  for  the  new 
best  fit  period  are  shown  in  Fig.  2.  In  the  following  we  use  the 
ephemeris  and  period  newly  derived  here. 

3.3.  Changes  in  the  spots 

Provided  the  newly  derived  ephemeris  represents  the  correct 
period  for  the  surface  structure  on  V410Tau,  the  minima  of 
phase  folded  lightcurves  for  each  observing  season  should  now 
coincide  with  0  =  0.  This  is  verified  in  Fig.  3  where  the 
phase  plots  of  the  yearly  averages  of  the  V  band  lightcurve 
of  V410Tau  since  1981  are  shown.  Most  of  the  photometry 
from  earlier  years  is  available  only  in  graphical  form,  and  the 
time  information  can  not  be  extracted  accurately  enough  for 
studies  of  the  1.87  d  periodicity  (original  references  Romano 
1975;  Rossiger  1981). 


3  The  T  Tauri  photometry  data  base  is  compiled  and  maintained  by 
W.  Herbst  and  available  at 
http : //www . astro . wesleyan . edu/~bi 11/ 
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Year  -  1900 


N 

Fig.  2.  O-C  diagram  for  the  seasonally  averaged  times  of  minimum  T0 
in  the  V  band  lightcurve  of  V410Tau  versus  rotational  cycle  num¬ 
ber  N.  Open  circles  -  for  the  period  given  by  P94,  and  filled  circles  - 
after  adapting  the  period  to  minimize  the  residuals.  Data  obtained  in 
Nov.  2001  were  used  as  reference  point. 


The  minima  of  all  years  included  in  our  determination  of 
the  period,  i.e.  1990  to  2001,  are  aligned  at  the  same  phase, 
as  expected.  The  minima  of  the  years  1986  and  1989  do  show 
a  shift.  This  phase  migration  seems  to  be  irregular  and  may 
indeed  indicate  changes  in  the  spot  location  or  distribution. 
Between  1981  and  1985  the  lightcurve  of  V410Tau  was  char¬ 
acterized  by  a  completely  different,  double-peaked  structure. 
Data  from  these  years  has  been  extensively  discussed  by  Herbst 
(1989)  who  showed  that  a  two-spot  model  provides  a  good  de¬ 
scription  of  the  visible  lightcurve. 

Changes  in  the  shape  of  the  lightcurve  reflect  variations 
of  the  structure  of  active  regions.  In  order  to  quantify  these 
changes  we  plot  the  time-evolution  of  amplitude,  mean,  min¬ 
imum,  and  maximum  of  the  V  band  lightcurve  in  Fig.  4.  Over 
the  years  the  lightcurve  of  V410  Tau  has  systematically  become 
more  stable  and  its  variability  more  regular. 

From  a  simple  harmonic  fit  to  the  mean  V  band  magnitude 
for  data  obtained  after  1990  (also  displayed  in  Fig.  4)  a  tenta¬ 
tive  cycle  length  of  5.4  yr  is  derived.  Inspection  of  the  photom¬ 
etry  in  other  bands,  also  available  at  the  T  Tauri  photometry 
data  base,  shows  the  same  5.4yr-pattern  in  the  B  band  with  a 
similar  amplitude,  supporting  the  idea  that  this  variation  can  be 
attributed  to  a  spot  activity  cycle.  The  data  in  the  R  and  I  band 
is  scarce  and  inhomogeneous  because  of  the  use  of  different 
filters  (Johnson  and  Cousin). 

However,  when  extrapolating  the  5.4yr-periodicity  to  ear¬ 
lier  years  (dotted  line  in  Fig.  4)  the  data  first  runs  out  of  phase 
with  the  5.4 yr  period  and  then  turns  into  completely  irregular 
behavior.  The  evolution  of  both  the  shape  of  the  lightcurve  and 
the  mean  brightness  are  suggestive  of  a  transition  towards  less 
active  behavior. 

4.  Simultaneous  X-ray  and  optical  data 

A  major  aim  of  this  campaign  was  to  study  simultaneous  X-ray 
and  optical  observations  in  order  to  look  for  correlations,  in 
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Fig.  3.  Seasonally  averaged  V  band  lightcurves  of  V410Tau 
from  1981  to  2001  phase  folded  with  the  new  ephemeris  given  in 
Table  2  (continued  on  the  next  page).  Data  are  extracted  from  the  data 
base  maintained  by  W.  Herbst.  The  data  for  2001  was  obtained  by  one 
of  us  (KG)  in  the  months  prior  to  our  coordinated  observing  campaign. 
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particular  those  related  to  the  spot  rotation  cycle.  To  this  end  the 
Chandra  observations  were  purposely  scheduled  to  cover  dif¬ 
ferent  rotational  phases  of  the  star.  If  X-ray  emission  is  related 
to  spotted  regions  the  maximum  of  the  X-ray  emission  should 
occur  at  the  minimum  of  the  optical  lightcurve,  i.e.  at  times 
when  the  spot  is  on  the  visible  hemisphere,  and  vice  versa. 
Therefore,  to  maximize  the  amplitude  of  the  expected  X-ray 
variability  we  observed  near  optical  maximum  and  minimum. 
In  addition  our  spectroscopic  optical  monitoring  provided  in¬ 
formation  on  the  time-evolution  of  the  Ha  line  and  the  radial 
velocity  (RV). 

The  V  band  lightcurve  obtained  during  our  monitoring  in 
November  2001  is  shown  in  the  lowest  panel  of  Fig.  5  phase- 
folded  with  the  new  period  and  ephemeris.  The  differential  pho¬ 
tometry  carried  out  in  the  Stromgren  system  was  transformed 
to  the  Johnson  system  with  help  of  the  absolute  Stromgren  pho¬ 
tometry  done  during  the  last  night.  We  point  out  that  the  pho¬ 
tometric  measurements  from  the  three  observing  sites  comple¬ 
ment  each  other  to  provide  nearly  full  coverage  of  the  rotational 
cycle  despite  the  relatively  short  duration  of  the  monitoring 
(11  d;  see  Fig.  1).  Photometric  data  in  the  UBRI  bands  was  also 
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Fig. 4.  Longterm  behavior  of  amplitude,  mean,  and  maximum  and  minimum  of  the  V  band  lightcurve  of  V410Tau  from  1978  to  2001.  Data 
points  represent  seasonal  averages.  The  numbers  on  the  top  i/-axis  stand  for  the  year  in  which  the  observing  season  started. 


obtained,  but  does  not  provide  new  information  on  the  phas¬ 
ing  of  the  1.87  d  cycle.  We  will  use  these  latter  lightcurves  for 
the  discussion  of  flares  in  an  accompanying  paper  (Fernandez 
et  al.,  in  prep.).  Some  of  these  flares  can  be  seen  in  the  V  band 
lightcurve  displayed  in  Fig.  5.  However,  our  aim  in  this  section 
is  to  examine  activity  parameters  for  a  possible  relation  to  the 
1 .87  d  cycle,  where  short-term  random  processes  such  as  flares 
are  not  of  interest. 

4.1.  X-rays 

The  mean  ACIS  count  rates  listed  in  Table  1  indicate  that  the 
X-ray  emission  of  V410  Tau  was  variable.  To  check  the  relation 
with  the  optical  photometry  we  generated  X-ray  lightcurves, 
and  phase-folded  them  with  the  1.871970d  period  derived 
in  Sect.  3.2.  The  result  for  the  two  observations  obtained  in 
November  2001  is  shown  in  Fig.  5.  A  rapid  rise  of  the  X-ray 
count  rate  took  place  near  optical  minimum  in  the  observation 
from  November,  probably  indicating  the  onset  of  a  flare.  This 
event  was  not  accompanied  by  simultaneous  observations  in 
the  optical,  but  just  a  few  hours  later  two  flares  were  observed 


from  Mt.  Maidanak,  hinting  at  a  possible  relation  between  the 
optical  and  X-ray  emission  sites. 

A  third  Chandra  pointing  was  obtained  in  March  2002 
near  optical  minimum,  i.e.  coinciding  in  phase  with  the  first 
November  observation.  As  this  observation  is  not  simultane¬ 
ous  with  the  optical  lightcurve  we  do  not  display  it  in  Fig.  5. 
The  March  data  does  not  show  significant  time-variability, 
and  its  count  rate  is  somewhat  higher  than  the  lowest  (quies¬ 
cent)  level  measured  in  the  November  observation  taken  at  the 
same  rotational  phase.  Thus,  while  the  two  November  pointings 
seem  to  suggest  that  the  X-ray  emission  depends  on  rotational 
phase,  this  conclusion  is  not  supported  by  the  observation  from 
March. 

4.2.  Ha 

The  Ha  line  is  in  the  spectral  range  of  both  the  intermediate- 
and  the  high-resolution  spectra  obtained  within  this  campaign. 
Ha  is  always  weak  in  emission.  We  display  the  equivalent 
width  Wna  measured  at  different  phases  of  the  photometric  ro¬ 
tation  cycle  in  the  second  panel  of  Fig.  5.  On  Nov.  20  both 
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Fig.  5.  Multi-wavelength  phase  plot  for  the  1.87  d  cycle  of  V410Tau.  From  top  to  bottom:  Radial  velocity,  equivalent  width  of  Ha,  X-ray  count 
rate,  and  V  band  photometry.  Except  for  some  data  points  in  the  radial  velocity  curve  all  measurements  are  obtained  quasi-simultaneously 
in  November  2001.  The  large  flare  visible  in  the  V  band  lightcurve  was  also  observed  spectroscopically,  but  the  equivalent  widths  of  Ha  are 
outside  the  range  displayed  in  this  figure  (see  Fernandez  et  al.,  in  prep,  for  a  detailed  discussion).  The  fast  increase  of  the  X-ray  count  rate  near 
0  =  0  (marked  with  an  arrow)  probably  indicates  the  beginning  of  a  flare  (see  Sect.  5). 


intermediate-  and  high-resolution  spectra  were  obtained,  and 
the  measured  equivalent  widths  are  in  good  agreement  within 
the  uncertainties.  The  error  bars  for  the  high-resolution  spectra 
are  based  on  the  assumption  of  an  accuracy  in  the  continuum 
normalization  of  10%.  The  typical  error  bar  for  the  values  from 
the  low-resolution  spectrum  are  -30%. 

On  Nov.  24  a  large  flare  erupted  which  we  monitored  simul¬ 
taneously  in  spectroscopy  and  photometry.  Note,  that  the  event 
is  not  seen  in  the  equivalent  width  curve  shown  in  Fig.  5  be¬ 
cause  the  values  for  Wna  measured  during  the  flare  are  higher 
than  the  displayed  plot  range. 


Due  to  poor  weather  conditions,  the  phase  coverage  of  the 
Her  data  is  limited.  While  the  equivalent  width  seems  to  un¬ 
dergo  (short-term)  changes  within  each  night,  no  clear  trend 
related  to  the  rotational  phase  can  be  found. 

4.3.  Radial  velocity  and  v  sin  i 

We  derived  the  RV  and  the  projected  rotational  veloc¬ 
ity  (v  sin  i)  from  the  optical  high-resolution  spectra  using  a 
cross-correlation  analysis  with  several  template  stars  of  spec¬ 
tral  class  K2,  K4  and  K5. 
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The  value  of  i>sint  is  obtained  by  comparing  the  spectra 
of  V410Tau  within  several  spectral  orders  to  that  of  a  broad¬ 
ened  template.  Amongst  the  templates  at  our  disposal,  the  K2  V 
star  HD  166620  broadened  to  a  velocity  of  74  +  3  k  m  s  1  gave 
the  best  fit.  This  broadening  agrees  with  the  estimates  of 
Vogel  &  Kuhi  (1981)  (76  +  lOkms-1),  Hartmann  et  al.  (1986) 
(70.9  +  9  km  s'1),  and  Hatzes  (1995)  (77  ±  1  kins"1). 

We  estimated  the  stellar  RV  from  the  data  acquired  in 
Nov.  2001.  This  was  done  by  averaging  measurements  that 
were  obtained  at  phases  0.5  <  (p  <  0.6.  The  four  spectra  ob¬ 
tained  in  this  phase  interval  correspond  to  the  maximum  in  the 
lightcurve,  and  therefore  their  shape  should  be  less  distorted 
by  the  spots.  Using  this  procedure  we  found  an  average  RV 
of  17.9  +  1.8kms_1  similar  to  the  stellar  RV  of  V410Tau  given 
in  the  literature  (1 8  l< m  s  ~ 1 ,  Herbig  &  Bell  1988). 

The  relative  shift  between  photospheric  lines  was  then 
measured  by  cross-correlating  the  spectra  of  V410Tau  with 
the  spectrum  of  the  template  HD  166620.  Three  different  or¬ 
ders  have  been  used,  with  wavelengths  centered  at  5400,  6100 
and  6400  A,  and  a  wavelength  coverage  of  80  A  in  each  order. 
A  parabolic  function  was  used  to  find  the  center  and  width  of 
the  correlation  peak,  and  the  errors  were  computed  from  the  fit¬ 
ted  peak  height  and  the  asymmetric  noise  as  described  by  Tonry 
&  Davis  (1979).  The  final  RV  are  an  average  of  those  obtained 
with  the  three  orders,  and  the  accuracies  are  the  combined  error 
bars  for  each  individual  measurement. 

We  performed  this  procedure  for  each  of  the  12  individ¬ 
ual  high-resolution  spectra  of  V410Tau  obtained  during  our 
campaign.  The  resulting  RV  curve  is  shown  in  the  top  panel 
of  Fig.  5  as  a  function  of  rotational  phase,  combined  with  a 
number  of  measurements  obtained  in  1993  and  presented  by 
Fernandez  &  Miranda  (1998)  added  here  to  improve  the  phase 
coverage. 

From  this  latter  data  set  we  used  only  the  spectrum  in 
the  range  6680-6730  A.  This  spectral  range  is  useful  because 
strong  photospheric  lines  such  as  the  Li  I  line  (6708  A)  are 
present.  In  a  first  step,  each  of  the  16  spectra  from  1993  was 
cross-correlated  with  the  first  spectrum  of  the  1993  series.  The 
position  of  the  maximum  of  the  cross-correlation  function  gives 
then  for  each  spectrum  the  average  shift  of  the  photospheric 
lines  relative  to  the  first  one.  To  be  able  to  combine  these  re¬ 
sults  with  the  RV  measurements  obtained  during  2001  we  took 
the  velocity  shift  at  rotational  phase  <p  =  0.5  as  a  zero-point.  At 
this  phase  the  photometry  shows  that  the  star  is  at  its  bright¬ 
est,  and  thus  the  spot(s)  or  a  large  part  of  it  is  not  visible. 
As  a  consequence,  we  expect  to  see  the  normal  stellar  photo¬ 
sphere,  and  the  velocity  shift  should  be  equal  to  the  stellar  RV, 
i.e.  AVrad  =  0. 

From  Fig.  5  it  is  readily  seen  that  the  RV  varies  from  w— 10 
to  + 1 0  km  s  1 ,  and  the  phase  shift  is  ~0.25  when  compared 
with  the  photometric  lightcurve.  This  phase  shift  can  be  un¬ 
derstood  if  the  cool  spot  is  responsible  for  the  distortion  of  ab¬ 
sorption  lines  (Vogt  et  al.  1987).  In  this  scenario,  the  RV  is 
large  when  the  spot  is  near  the  limb  and  close  to  zero  when 
the  spot  is  face-on  (0  «  1.0)  or  occulted  by  the  star  (</>  »  0.5). 
This  result  shows  how  magnetic  stellar  activity  can  affect  the 
spectroscopic  search  for  very  low-mass  companions  orbiting 
around  PMS  stars. 


5.  X-ray  properties 


We  extracted  an  X-ray  spectrum  for  each  of  the  three  Chandra 
pointings.  The  observation  from  Nov.  16/17  shows  a  pro¬ 
nounced  increase  in  count  rate  at  the  end  of  the  exposure 
(discussed  above  in  Sect.  4.1).  If  this  behavior  in  the  X-ray 
lightcurve  is  due  to  a  flare,  it  is  expected  that  the  spectrum 
changes  as  a  consequence  of  coronal  heating.  To  unveil  any 
eventual  spectral  variability  in  the  data  we  split  the  Nov.  16/17 
observation  into  a  quiescent  and  flaring  part,  separating  the  data 
at  JD  2  452  230.4870  (marked  with  an  arrow  in  Fig.  5).  No  sys¬ 
tematic  trend  is  seen  in  the  temporal  behavior  of  the  other  two 
Chandra  pointings.  Therefore,  we  think  of  both  as  representing 
the  quiescent,  i.e.  non-flaring,  state  of  V410Tau. 

As  outlined  in  Sect.  2.3  we  avoided  excessive  photon  pile- 
up  by  placing  V410Tau  at  an  off-axis  angle  of  ~4'.  To  check 
whether  any  remaining  pile-up  affects  the  spectrum  the  source 
extraction  area  was  varied:  next  to  the  spectrum  from  a  4"  ra¬ 
dius  centered  on  the  wavdetect  position  of  V410Tau  we  ex¬ 
tracted  the  spectrum  from  an  annulus  where  we  excluded  the 
inner  source  region.  In  the  outer  portions  of  the  PSF  the  count 
rate  is  smaller  and  no  pile-up  is  expected.  Visual  inspection 
does  not  show  any  distortion  of  the  full-source  spectrum  with 
respect  to  the  one  extracted  from  a  (with  certainty  pile-up  free) 
annulus.  The  supposition  that  pile-up  is  negligible  was  con¬ 
firmed  when  we  fitted  both  spectra  with  the  same  model  and 
found  no  significant  differences  in  the  spectral  parameters.  In 
the  following  we  describe  the  details  about  our  spectral  model. 

Spectral  fitting  was  done  in  the  XSPEC  environment  (ver¬ 
sion  11.2.0).  To  make  up  for  the  continuous  degradation  of 
the  ACIS  quantum  efficiency  we  applied  the  acisabs  model 
to  the  auxiliary  response  file  (arf)  before  loading  the  data 
into  XSPEC. 

We  applied  the  MEKAL  model  (Mewe  et  al.  1985)  for  ther¬ 
mal  emission  from  an  optically  thin,  hot  plasma  plus  a  photo¬ 
absorption  term.  Comparison  of  the  three  quiescent  spectra 
(first  part  of  Nov.  16/17,  Nov.  19/20,  and  Mar.  7)  shows  that 
their  spectral  shape  is  very  similar.  Separate  fitting  led  to  in¬ 
distinguishable  parameters.  Therefore,  we  present  here  the  re¬ 
sult  of  a  joint  modelling  of  all  three  quiescent  spectra  which 
improves  the  statistics.  To  represent  the  observed  small  offsets 
in  the  flux  of  the  three  spectra  we  allowed  for  an  independent 
normalization  constant. 

An  acceptable  fit  (yjL  ~  1  and  flat  residual)  requires  a  min¬ 
imum  of  three  MEKAL  components,  and  free  elemental  abun¬ 
dances.  The  individual  spectra,  the  best  fit  model  and  residuals 
are  displayed  in  Fig.  6.  In  Table  3  we  summarize  the  spectral 
parameters  observed  during  the  quiescent  state. 

To  examine  whether  the  rise  in  the  lightcurve  on  Nov.  16/17 
represents  a  heating  event  we  make  use  of  hardness  ratios. 
Hardness  ratios  are  defined  as  follows: 


HR  = 


B\  -B2 
BI+B2 


(1) 


where  B\  and  B2  denote  the  count  rates  in  a  hard  and  a  soft 
band,  respectively.  We  defined  two  hardness  ratios  for  the  spec¬ 
tral  range  of  ACIS:  HR  1  based  on  the  range  from  0. 2-1.0  keV 
as  the  soft  band,  and  1.0-2.0keV  as  the  hard  band.  And  HR2 
with  1.0-2.0keV  representing  the  soft  band  and  2.0-8.0keV 
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Fig.  6.  Quiescent  X-ray  spectrum  of  V410Tau  during  three  Chandra 
exposures.  A  flare  has  been  eliminated  from  Obs.  Seq.#  200130.  The 
data  are  overlaid  by  the  best  fit  model  from  Table  3.  The  residuals  for 
each  observation  are  shown  in  separate  panels. 


representing  the  hard  band.  In  Fig.  7  these  hardness  ratios  are 
plotted  for  the  four  time  segments  introduced  above.  A  harden¬ 
ing  of  the  spectrum  with  a  significance  of  ~2  cr  is  apparent  dur¬ 
ing  the  second  part  of  the  Nov.  16/17  observation.  This  may  be 
seen  as  an  indication  that  indeed  V410Tau  underwent  a  flare, 
because  flares  are  a  result  of  (magnetic)  heating. 


6.  Discussion 

We  have  carried  out  an  intensive  coordinated  monitoring  cam¬ 
paign  in  the  optical  and  X-ray  wavelength  ranges  with  the 
aim  to  study  correlations  between  the  photometric  rotation  cy¬ 
cle  of  V410Tau  and  different  activity  diagnostics.  Combining 
our  new  data  with  historic  photometric  measurements  we  re¬ 
examined  the  variability  in  the  optical  lightcurve  of  V410Tau 
on  various  timescales:  (i)  spot  (rotation)  cycle,  and  (ii)  long¬ 
term  (activity)  cycle. 

6.1.  Variability  in  the  spot  rotation  cycle 

Optical  photometric  observations  were  performed  at  three  sites 
around  the  globe,  thus  providing  complete  phase  coverage  of 
the  1.87  d  spot  cycle  despite  the  short  monitoring  time  of  1 1  d. 
This  has  allowed  us  to  measure  the  time  of  the  minimum  in  the 
V  band  lightcurve  of  V410  Tau  with  high  precision.  Combining 
this  measurement  with  lightcurves  from  the  years  1990  to  2001 
we  derived  an  update  of  the  rotational  period.  The  need  for  a 
revision  of  the  period  was  indicated  by  a  systematic  shift  of  the 


Table  3.  X-ray  spectral  parameters  of  V4 10  Tau  derived  from  a  joint  fit 
to  the  three  ACIS  observations  that  represent  the  stars’  quiescent  state. 
A  constant  normalization  factor  was  applied  to  allow  for  the  observed 
difference  in  X-ray  brightness  between  the  three  data  sets:  (V200190  =  1, 
N200191  =  0.91  ±0.04,  N200130  =  0.74±0.04.  Errors  are  90%  confidence 
levels. 


Parameter  Unit 

Fit  result 

Nh 

[1021  cnr3] 

0  97+0-37 
-0.32 

kTx 

[keV] 

0  24+0  03 
-0.02 

kT2 

[keV] 

0  Q9+0.09 
yj‘yj-0.07 

kT3 

[keV] 

9  1  O+0.43 
Z-i8-0.28 

EM[ 

[1053  cnr3] 

3  02+101 

em2 

[1053  cnr3] 

9  99+1.26 
-0.80 

EM3 

[1053  cnr3] 

9  1  q+0.45 

1  -0.77 

Z 

[Zo] 

0  20+008 
-0.06 

Wed  (d-°-f0 

1.30  (330) 

HR  1 


Fig.  7.  ACIS  hardness  ratios  for  V410Tau.  Flare  and  quiescent  state 
for  observation  200130  are  labeled  with  “F”  and  “Q”,  respectively. 
Error  bars  denote  lcr  uncertainties. 


times  of  minimum  in  seasonally  averaged  optical  lightcurves 
folded  with  the  most  recent  ephemeris  for  V410Tau  given 
by  P94.  We  attribute  this  to  the  fact  that  P94  based  their  pe¬ 
riod  determination  on  data  acquired  between  1986  and  1992, 
where  the  time  of  minimum  seems  to  have  varied  erratically. 
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The  new  value  for  the  period  is  slightly  smaller  than  the  period 
given  by  P94. 

An  alternative  explanation  for  the  monotonic  ally  increasing 
shift  of  the  time  of  minimum  light  observed  over  the  last  10  yrs 
is  latitudinal  migration  of  spots  on  the  differentially  rotating 
surface.  However,  making  use  of  the  differential  rotation  pa¬ 
rameter  of  V410  Tau  known  from  Doppler  imaging  we  find  that 
the  associated  change  in  period  would  correspond  to  a  quite 
large  movement  of  the  spot  on  the  star,  which  is  not  supported 
by  the  photometry  and  Doppler  images.  Therefore  we  consider 
this  scenario  unlikely  as  an  explanation  for  the  systematic  ef¬ 
fect  observed  since  1990.  Nevertheless,  spot  migration  may  be 
responsible  for  the  irregular  shifts  of  the  time  of  minimum  ob¬ 
served  in  the  years  before.  Unfortunately,  no  Doppler  images 
are  available  for  these  early  years.  But  spot  models  adapted 
to  the  photometric  data  imply  indeed  changes  in  the  surface 
distribution  (latitude,  longitude,  and  filling  factor)  of  the  spots 
(Herbst  1989;  Bouvier  &  Bertout  1989). 

In  our  series  of  high-resolution  spectra  we  detected 
RV  variations  which  can  be  explained  by  the  distortions  that 
a  spot  induces  onto  the  absorption  line  profiles  of  a  rotating 
star.  In  our  RV  curve  we  included  data  obtained  in  1993  next  to 
our  recent  monitoring  from  Nov.  2001.  The  RV  measurements 
of  these  two  years  agree  very  well  in  both  phase  (if  folded  with 
the  new  period)  and  amplitude  demonstrating  that  the  spot  dis¬ 
tribution  has  not  changed  significantly  over  the  last  decade.  The 
half-amplitude  of  the  measured  variability  is  AVmd  ~  9  km  s'. 
The  RV  curve  of  V410Tau  was  also  measured  by  Welty  & 
Ramsey  (1995).  Their  data  showed  an  indication  for  an  increase 
in  amplitude  from  ~6  to  ~8  km  s  1  from  1992  to  1993,  while 
the  overall  shape  remained  stable.  These  small  amplitude  varia¬ 
tions  are  consistent  with  the  changes  in  the  V  band  photometry 
of  the  same  years,  and  might  be  related  to  longterm  changes  in 
the  spot  size. 

In  recent  years  detailed  RV  studies  with  high-precision 
have  been  carried  out  for  mostly  solar-type  stars  with  the  aim  of 
detecting  low-mass  (planetary)  companions  around  them  (see 
e.g.  Queloz  et  al.  2001  for  a  review).  However,  in  active  stars 
spot-related  variability  may  dominate  the  RV  time  series,  such 
that  the  detection  of  planets  is  impeded.  Therefore,  it  is  im¬ 
portant  to  assess  the  amount  of  perturbation  induced  by  star 
spots.  Saar  &  Donahue  (1997)  calculated  the  maximum  pertur¬ 
bation  of  the  RV  that  a  spot  distribution  with  a  filling  factor  /s 
induces  onto  a  star:  As  «  6.5  X  jf-9  x  v  sin  i.  By  their  defini¬ 
tion  /s  characterizes  the  inhomogeneous  part  of  the  spot  pat¬ 
tern.  If  we  apply  their  Eq.  (1)  (given  above)  to  V410Tau  we 
derive  /s  ~  29%,  an  almost  perfect  match  to  the  actual  size  of 
the  spot  on  V410Tau  inferred  from  Doppler  imaging  (~30%; 
Joncour  et  al.  1994).  Our  estimate  represents  the  first  attempt 
to  extend  the  Saar-Donahue  relation  to  fast  rotating  (and  very 
active)  stars,  and  we  conclude  that  this  relation  seems  to  hold 
also  in  this  regime.  Furthermore,  our  result  implies  that  the  spot 
on  V410Tau  is  significantly  non-uniform.  Indeed,  Doppler  im¬ 
ages  indicate  that  the  major  spot  is  near  the  pole  but  not  quite 
identical  with  a  “polar  cap”. 

There  seems  to  be  no  other  component  in  the  RV  data  be¬ 
sides  the  spot  induced  variability.  We  use  the  uncertainties  in 
the  individual  RV  measurements  (~3kms  *)  to  determine  an 


upper  limit  for  the  mass  of  a  short-period  binary  companion 
possibly  hidden  in  the  data.  Assuming  a  circular  orbit  for  the 
hypothesized  companion,  we  can  exclude  such  objects  with  a 
mass  of  m  sin  i  >  0.03  MQ  and  a  period  of  10  d  or  less.  Due 
to  the  short  timescale  of  our  monitoring  no  information  can  be 
obtained  about  longer  period  companions. 

6.2.  Detection  of  an  activity  cycle? 

Inspection  of  the  V  band  lightcurves  of  V410Tau  observed 
during  the  last  two  decades  indicates  that  its  behavior  has  be¬ 
come  systematically  more  stable  and  its  variability  more  reg¬ 
ular:  The  first  photometric  data  set  which  provides  precise  in¬ 
formation  on  both  the  observing  time  and  the  magnitude,  with 
good  coverage  in  phase,  is  from  1981.  Around  this  time  a  dis¬ 
tinct  double-peaked  structure  evolved,  which  smoothly  trans¬ 
formed  itself  into  a  single-peaked  lightcurve  over  the  follow¬ 
ing  five  years.  This  change  went  along  with  an  increase  in  the 
amplitude  by  a  factor  of  2-3.  Since  then  the  simple  near-to 
sinusoidal  shape  of  the  lightcurve  has  persisted,  only  subject 
to  comparatively  small  variations.  The  most  obvious  variation 
within  the  last  ten  years  is  the  systematic  pattern  seen  in  the 
mean  amplitude.  From  a  simple  harmonic  fit  to  the  seasonal 
averages  of  the  V  band  magnitude  of  the  years  1990  to  2001  a 
periodic  variation  of  5.4 yr  length  is  suggested  that  might  rep¬ 
resent  an  activity  cycle. 

Baliunas  et  al.  (1996)  have  outlined  a  connection  between 
the  ratio  of  cycle  and  rotation  periods,  Pcyc/Pmt,  with  the  dy¬ 
namo  number,  D ,  measuring  the  efficiency  of  the  stellar  dy¬ 
namo.  According  to  stellar  dynamo  theory  D  is  directly  pro¬ 
portional  to  l/Prot-  From  observations  of  chromospherically 
active,  slowly  rotating  stars  they  found  that  the  empirical 
slope  in  the  lg  (PCyc/Prot)  versus  lg(l/Prot)  diagram  is  =s0.74. 
Interpreting  the  5 .4  yr- variability  of  V4 1 0  Tau  as  a  cycle  period 
we  find  that  its  location  in  the  lg  (Pcyc/Pmt)  ~  lg  (1/Prot)-  dia¬ 
gram  is  consistent  with  the  extrapolation  of  the  line  identified 
by  Baliunas  et  al.  (1996)  into  the  regime  of  fast  rotating  stars. 
This  result  is  remarkable  because  it  seems  to  indicate  that  the 
dynamo  process  on  PMS  stars  is  very  similar  to  that  on  more 
evolved  stars.  To  the  best  of  our  knowledge  V410Tau  is  the 
first  PMS  star  investigated  in  this  respect. 

To  extend  the  comparison  with  the  evolved  stars  we  exam¬ 
ined  the  position  of  V410Tau  in  the  lg  («cyc/£2rot)  -lgR0'  -  di¬ 
agram,  where  Rq  =  2rCOnv^rot  is  the  Rossby-number,  and  Tconv 
the  convective  turnover  time.  Saar  &  Brandenburg  (1999)  de¬ 
rived  TConv  for  their  sample  of  chromospherically  active  stars 
(mostly  objects  from  the  Mt.  Wilson  Call  H+K  survey)  from 
the  models  of  Gunn  et  al.  (1998),  and  identified  three  branches 
in  this  plot:  inactive,  active,  and  superactive  stars.  To  esti¬ 
mate  TConv  for  V410Tau  we  compare  its  stellar  parameters 
to  the  PMS  calculations  by  Ventura  et  al.  (1998),  and  find 
lg  Tconv  Id]  ~  2.3.  Because  we  have  found  agreement  within 
0.2  dex  between  zero-age  main  sequence  convective  turnover 
times  estimated  from  models  by  Ventura  et  al.  (1998)  and  those 
from  models  by  Gunn  et  al.  (1998),  we  can  now  make  valid 
comparisons  between  V410Tau  and  main-sequence  stars.  We 
find  that  the  location  of  V410Tau  is  not  consistent  with  any  of 
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the  regions  defined  by  the  evolved  stars  in  the  lg  («cyc/£2rot)  - 
Ig  R{.  1  -  diagram.  This  is  presumably  due  to  its  fast  rotation 
combined  with  a  large  convective  turnover  time,  which  is  a 
consequence  of  its  PMS  nature.  Indeed,  most  PMS  stars  -  typ¬ 
ically  characterized  by  PI0t  <  lOd  and  rCOnv  ~  60...400d  - 
are  expected  to  lie  to  the  right  of  the  active/inactive  branches 
in  this  diagram  (lgRl1  >  2),  and  below  these  branches.  Given 
their  generally  fast  rotation  only  very  short  cycle  periods  would 
place  them  at  the  extension  of  these  lines. 

V410Tau  is  younger  (~1  Myr)  than  all  other  stars  for  which 
activity  cycles  have  been  reported  so  far,  representing  a  unique 
test  case  for  dynamo  action  on  the  PMS.  Among  the  stars  that 
seem  to  display  cyclic  behavior  it  comes  closest  to  the  young 
solar-analogs,  such  as  AB  Dor,  EK  Dra,  and  LQ  Hya,  which 
are  single  stars  characterized  by  an  age  of  ~50-80Myr,  spec¬ 
tral  type  of  early  K,  and  fast  rotation  ( PI0t  <  3d).  A  cycle 
length  for  AB  Dor  of  5.3  yr  was  given  by  Amado  et  al.  (2001). 
Berdyugina  et  al.  (2002)  identified  three  cycles  of  5.2,  7.7, 
and  15  yr  duration  in  LQHya.  The  photometry  of  EKDra  in¬ 
dicates  longterm  fading  over  the  last  35  years  (Frohlich  et  al. 
2002),  while  Saar  &  Brandenburg  (1999)  have  predicted  cycle 
periods  of  1 .4  and  39  yr  for  this  star. 

The  observation  of  activity  cycles  in  young  stars  is  of 
paramount  importance  for  the  understanding  of  the  dynamo  op¬ 
erating  in  these  objects.  The  solar-type  crQ-dynamo  is  thought 
to  be  localized  in  the  overshoot  layer  at  the  bottom  of  the  con¬ 
vection  zone.  But  in  stars  with  deep  convective  envelope  a  “dis¬ 
tributed”  dynamo  located  throughout  the  convective  layer  may 
take  over.  Furthermore,  both  observations  and  theory  agree  in 
that  differential  rotation  is  suppressed  in  fast  rotating  young 
stars  (Henry  et  al.  1995;  Ktiker  &  Rudiger  1997).  Therefore, 
the  o'-effect  should  dominate  over  the  rotational  shear,  and  a 
pure  G-O-dynamo  is  unlikely  to  hold.  The  type  and  stability  of 
the  dynamo  solutions  depend  critically  on  the  strength  of  differ¬ 
ential  rotation.  For  small  values  of  the  differential  rotation  non- 
axisymmetric  modes  are  preferred  (Moss  et  al.  1995).  These 
modes  seem  not  to  oscillate  (Ktiker  &  Rudiger  1999)  and  are 
difficult  to  reconcile  with  the  possible  observation  of  an  activity 
cycle  on  V410Tau.  In  fact,  the  apparent  absence  of  such  cycles 
in  PMS  stars  is  usually  taken  as  evidence  for  the  action  of  a 
non-solar  dynamo. 

However,  Kitchatinov  et  al.  (2001)  have  shown  that  a  tran¬ 
sition  of  the  dynamo  mode  takes  place  at  an  age  of  ~5-10  Myr, 
such  that  in  the  older  stars  an  axisymmetric  oscillating  field  is 
preferred.  The  age  of  V410Tau  is  close  to  this  critical  range, 
indicating  that  it  may  be  a  transition  object.  This  might  explain 
the  presence  of  an  activity  cycle  (representative  for  oscillating 
fields)  in  conjunction  with  long-lived  spots  (indicating  very  sta¬ 
ble  field  structures). 

Our  present  knowledge  draws  a  complicated  picture  for  the 
magnetic  activity  of  V410Tau,  and  PMS  stars  in  general.  We 
stress  that  our  conclusions  rely  on  the  observation  of  just  one 
tentative  cycle  period.  The  systematic  pattern  in  the  season¬ 
ally  averaged  mean  magnitude  seems  to  have  its  onset  a  few 
years  before  1990,  at  the  same  time  when  the  single-peaked 
lightcurve  developed.  Earlier  than  — 1986  the  lightcurve  shows 
no  signs  for  cyclic  variability.  Indeed,  in  years  of  low  amplitude 
the  V  band  magnitude  assumes  an  intermediate  level,  while  a 


value  near  maximum  would  be  expected  if  a  decrease  in  spot 
size  and  number  during  a  cycle  minimum  were  responsible. 

6.3.  Spot  rotation  cycle  and  activity  diagnostics 

Simultaneously  with  the  optical  photometric  lightcurve  we 
have  acquired  Chandra  X-ray  observations  and  optical  spec¬ 
troscopy.  Chandra  has  targeted  V410Tau  twice  during  our 
campaign,  at  minimum  and  maximum  optical  brightness,  re¬ 
spectively.  Different  count  levels  were  found,  but  their  relation 
to  the  rotation  cycle  was  not  confirmed  by  a  following  measure¬ 
ment  carried  out  some  months  later  that  showed  an  intermedi¬ 
ate  count  rate  although  obtained  at  the  same  phase  as  one  of  the 
observations  from  Nov.  2001.  A  seeming  lack  of  a  correlation 
between  X-ray  and  optical  emission  was  already  pointed  out 
by  our  analysis  of  archived  ROSAT  data  (Stelzer  et  al.  2002): 
we  found  that  the  count  rates  varied  from  one  observation  to 
the  other,  however,  without  clear  relation  to  the  rotation  cycle. 
Similarly  the  Her  equivalent  width  seems  not  to  show  a  trend 
related  to  the  rotational  phase.  Unfortunately,  we  were  not  able 
to  obtain  full  phase  coverage  in  the  optical  spectroscopy  due  to 
poor  weather  conditions. 

Reports  on  coordinated  multi-wavelength  monitoring  of 
PMS  stars  are  scarce  in  the  literature.  Multi-wavelength  obser¬ 
vations  of  TTS  in  the  Taurus  star  forming  region  using  ROSAT 
jointly  with  optical  telescopes  have  revealed  flares  in  Balmer 
lines  and  in  X-rays  (Guenther  et  al.  2000).  But  due  to  unfor¬ 
tunate  conditions  none  of  these  events  was  observed  simulta¬ 
neously  in  both  the  optical  and  the  X-ray  range.  In  the  same 
study  a  weak  correlation  between  the  X-ray  and  Her  emission 
was  seen  for  the  wTTS  V773  Tau,  suggestive  of  a  relation  be¬ 
tween  the  emission  sites.  Simultaneous  X-ray  and  optical  ob¬ 
servations  of  BP  Tau  were  discussed  by  Gullbring  et  al.  (1997). 
No  signs  for  any  correlation  between  the  optical  and  the  X-ray 
emission  was  seen:  Two  optical  flares  had  no  counterpart  in 
X-rays,  and  the  data  set  did  not  allow  them  to  examine  varia¬ 
tions  related  to  the  rotation  of  the  star.  However,  BP  Tau  is  a 
cTTS,  such  that  variability  (in  both  optical  and  X-rays)  may  be 
induced  by  accretion.  Therefore,  it  may  not  be  directly  com¬ 
parable  to  V410Tau  which  is  a  non-accreting  wTTS  where 
all  variability  should  be  related  to  magnetic  activity  similar  to 
more  evolved  stars. 

Recently,  Jardine  et  al.  (2002)  have  modeled  the  X-ray 
emission  of  ABDor,  for  which  Kiirster  et  al.  (1997)  did  not 
find  any  evidence  for  rotational  modulation  of  the  X-ray  emis¬ 
sion  during  monitoring  with  ROSAT.  The  model  of  Jardine 
et  al.  (2002)  is  based  on  Zeeman  Doppler  maps  of  the  surface 
magnetic  field  structure,  and  predicts  little  rotational  modula¬ 
tion  of  the  X-ray  emission  due  to  the  extended  structure  and/or 
high  latitudes  of  coronal  features.  AB  Dor  is  characterized  by 
dark  spots  at  all  latitudes  (Donati  &  Collier  Cameron  1997). 
On  V410Tau  the  dominant  spot  seems  to  be  located  near  the 
pole,  although  not  completely  symmetrically  (Hatzes  1995). 
But  the  absence  of  rotational  modulation  in  the  X-ray  emis¬ 
sion  of  V410Tau  is  consistent  with  a  symmetric  distribution  of 
X-ray  flux  with  respect  to  the  rotation  axis.  Therefore,  the  spot 
may  not  be  the  main  site  of  X-ray  production. 
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Table  4.  X-ray  and  Ha  flux  and  respective  luminosities  of  V410Tau 
during  the  campaign  in  November  2001.  X-ray  emission  was  mea¬ 
sured  in  the  0.4-8  keV  energy  interval. 


Chandra 

Obs-ID 

/x 

[erg  cm-2  s~'] 

IgTx 

[erg/s] 

2001 3 0Q 

1.3  x  10-12 

30.5 

200190 

1.8  x  10-12 

30.6 

200191 

1.6  x  10-12 

30.5 

Rc 

/h  a 

lgl-Ha 

[A] 

[mag] 

[ergcnL2  s~*] 

[erg/s] 

Min 

1.15 

10.35 

1.9  x  10~13 

29.6 

Max 

2.45 

9.96 

5.7  x  10-13 

30.1 

6.4.  Relation  between  X-ray  and  Ha  flux 

Coronal  X-rays  and  chromospheric  Her  emission  are  thought  to 
be  produced  by  the  same  magnetic  heating  mechanisms  which 
are  thought  to  involve  either  magnetic  waves  (Goossens  1994) 
or  magnetic  reconnection  processes  (Priest  &  Forbes  2000)  as 
heating  agent.  Therefore  it  is  expected  to  find  a  correlation  be¬ 
tween  chromospheric  and  coronal  emission  when  comparing 
stars  at  different  activity  levels.  In  fact  Fleming  et  al.  (1988) 
and  Doyle  (1989)  have  found  such  a  correlation  for  a  sam¬ 
ple  of  dMe  flare  stars.  For  the  PMS  no  comparable  studies 
exist  mainly  because  of  a  lack  of  Ha  flux  measurements.  As 
we  will  show  in  an  accompanying  paper  (Fernandez  et  al.,  in 
prep.)  the  non-periodic  component  of  the  photometric  variabil¬ 
ity  of  V410Tau  is  reminiscent  of  flare  stars,  suggesting  that 
a  comparison  of  its  activity  to  that  latter  class  of  objects  is 
justified. 

We  estimated  the  X-ray  flux  of  V410Tau  from  the 
ACIS  spectrum.  The  result  is  given  in  Table  4.  Our  optical 
spectra  are  not  flux-calibrated.  But  the  Ha  flux  can  be  calcu¬ 
lated  from  Wna  with  help  of  the  photometry  that  we  can  use 
to  determine  the  continuum  flux,  because  of  the  very  regu¬ 
lar  pattern  of  the  lightcurve.  The  Wua  was  measured  on  the 
low-resolution  spectra,  taking  into  account  the  photospheric 
absorption.  We  multiplied  the  Wh„  by  the  specific  flux  of  the 
Rc  band  (in  erg  s  1  cm  2  A  1 )  and  used  the  Hipparcos  dis¬ 
tance  for  V410Tau  of  136 pc  (Wichmann  et  al.  1998)  to  com¬ 
pute  LHn.  The  specific  flux  for  a  star  with  Rc  =  Omag  can  be 
obtained  from  Rydgren  et  al.  (1984),  taking  into  account  the 
central  wavelength  of  the  band.  We  compute  the  Ha  flux  for  the 
minimum  and  maximum  equivalent  width  measured  in  the  low- 
resolution  spectra  during  the  quiescence  of  the  star.  The  corre¬ 
sponding  R  band  brightness  is  extracted  from  the  lightcurve  at 
the  same  rotational  phase.  The  minimum  and  maximum  quies¬ 
cent  Ha  flux  during  our  campaign  derived  in  this  way  is  tabu¬ 
lated  in  Table  4. 

For  the  M  dwarfs  studied  the  X-ray  luminosity  seems  to  be 
somewhat  higher  than  the  Ha  luminosity.  Fleming  et  al.  (1988) 
found  the  relation  LHn  ~  0.7 Lx.  Similarly,  Hawley  et  al.  (1996) 
quoted  values  between  ~0.2...  1  for  Ig  (LyJ for  both  field 
M  dwarfs  and  the  zero-age  main  sequence  cluster  IC  2602.  The 


values  we  derive  for  the  X-ray  and  Ha  luminosity  of  V410Tau 
are  in  good  agreement  with  these  relations,  suggesting  a  tight 
connection  between  the  activity  of  this  PMS  star  and  that  of 
MS  flare  stars. 


7.  Conclusion 

We  have  presented  the  results  from  a  coordinated  multi¬ 
wavelength  observing  campaign  for  the  wTTS  V410Tau  aim¬ 
ing  to  disentangle  the  role  that  the  various  atmospheric  layers 
play  in  magnetic  activity.  A  multi-wavelength  approach  is  es¬ 
sential  in  establishing  the  structure  and  relation  between  the 
emission  sites. 

Our  observations  of  V410Tau  can  be  summarized  as 
follows: 

-  An  update  of  the  period  and  ephemeris  of  the  1.87  d  cy¬ 
cle  representing  the  rotation  period  of  the  spot  dominating 
since  ~1990. 

-  The  non-detection  of  rotational  modulation  in  the  X-ray 
and  Ha  lightcurves  suggests  either  very  long  coronal  loops 
or  that  X-ray  emission  is  not  concentrated  in  the  major  spot, 
which  extends  over  the  pole  of  the  star  but  not  in  a  symmet¬ 
rical  way. 

-  The  first  tentative  detection  of  an  activity  cycle  on  a 
PMS  star  seems  to  indicate  that  the  field  generation  mecha¬ 
nism  that  drives  activity  in  young  stars  may  be  more  similar 
to  the  standard  solar-type  dynamo  than  widely  believed. 

-  The  detection  of  flares  in  both  the  optical  and  the  X-ray 
regime,  as  well  as  the  observed  relation  between  X-ray  and 
Ha  flux  put  V410  Tau  in  close  vicinity  to  the  more  evolved 
dMe  flare  stars. 

Whether  the  characteristics  of  V410Tau  are  peculiar  or  typical 
for  the  PMS,  or  whether  it  represents  a  rare  kind  of  a  transition 
object  remains  open  to  date.  We  would  like  to  encourage  sys¬ 
tematic  multi-wavelength  studies  of  other  young  stars  that  will 
eventually  put  the  nature  of  magnetic  activity  on  the  PMS  in  a 
clear  context  to  the  case  of  the  Sun  and  solar-like  stars. 
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